Analytic Kerr black hole lensing for equatorial observers in the strong deflection limit 
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In this paper we present an analytical treatment of gravitational lensing by Kerr black holes 
in the limit of very large deflection angles, restricting to observers in the equatorial plane. We 
accomplish our objective starting from the Schwarzschild black hole and adding corrections up to 
second order in the black hole spin. This is sufficient to provide a full description of all caustics and 
the inversion of lens mapping for sources near them. On the basis of these formulae we argue that 
relativistic images of Low Mass X-ray Binaries around Sgr A* are very likely to be seen by future 
X-ray interferometry missions. 
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I. INTRODUCTION 

General relativity predicts that light rays passing close 
to a black hole suffer gravitational lensing, so that an ob- 
server almost aligned with the line connecting a source 
and a black hole sees two images of the original source. 
These images are due to small deviations of photons that 
pass far enough from the black hole to allow a weak field 
approximation of the metric tensor. However, already 
Darwin in 1959 noticed that photons passing very close 
to a black hole may suffer much larger deflections with- 
out falling into the event horizon In principle, an 
observer situated on the line connecting the source and 
the black hole, besides the two classical weak field im- 
ages, would detect two infinite series of images very close 
the black hole. These images are produced by photons 
making one or more complete loops around the black hole 
before re-emerging in the observer direction. Of course, 
these relativistic images are largely demagnified w.r.t the 
original source and for some time they just remained a 
mathematical curiosity of General Relativity. 

Nevertheless, things changed after the great progress of 
interferometric techniques and the widely-accepted opin- 
ion that the radio-source Sgr A* in the Galactic center 
actually hosts a supermassive black hole of 3.61 x 10^ so- 
lar masses (for a review see Q). These facts motivated 
Virbhadra & Ellis to propose that this black hole may be 
an ideal candidate for the generation of relativistic im- 
ages of sources eventually passing behind it Q. In fact, 
the angular radius of the shadow of Sgr A* is predicted 
to be 23 /ias, which is comparable to the best resolution 
achieved in the millimeter band (18 /ias 5]). A complete 
imaging in the sub- mm band was suggested in Ref. 6J. 
Future space missions in the infrared and in the X-rays 
may reach even higher resolutions (for a complete discus- 



sion, see Ref. 0). 

A new field of Gravitational lensing was definitively 
opened and several authors proposed alternative meth- 
ods to overcome the evident difficulties of full general rel- 
ativity calculations of geodesies which typically result in 
cumbersome equations and heavy numerical integrations 
[^0. However, Darwin himself proposed a surprisingly 
easy formula for the positions of the relativistic images 
generated by a Schwarzschild black hole 1] . This formula 
and its consequences were later discussed or re-discovered 
several times [Hini IH El before Virbhadra & Ellis 
proposal. After that work, it was revived in Ref. [l3 |. 
where it was called the strong field limit of the deflection 
angle. It was then extended to Reissner-Nordstrom black 
holes in Ref. [isll and applied to microlensing by Sgr A* 
by Fetters [la|- In this paper, as suggested by Perhck 
|9| , we shall revise this terminology, referring to a Strong 
Deflection Limit (SDL), since an infinite deflection angle 
is not necessarily related to a large curvature. This can 
be realized by considering a very large black hole. The 
minimum distance reachable by a photon without being 
captured is of the same order of the horizon radius. The 
Riemann invariant RafijsR"'^^^ , evaluated in the inner 
region probed by the photon, scales as the curvature at 
the horizon, i.e. l/M'*. Increasing the mass of the black 
hole, the curvature felt by the photon becomes arbitrarily 
low, even if its defiection may be large. So, it is more cor- 
rect to speak of a strong deflection limit without referring 
to the curvature. 

The power of the SDL expansion of the deflection angle 
became evident when its universality was demonstrated 
in Ref. ^3- class of spherically symmetric black 

holes leads to the same SDL expansion; the coefficients 
of this expansion depend on the specific class of the black 
hole, representing a sort of identity card, from which all 
the parameters of the black hole can be extracted. By ob- 
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serving the relativistic images of a gravitational lensing 
event induced by a black hole, it is possible, in princi- 
ple, to deduce all its parameters and properties. Since 
this could also provide the key to discriminate between 
General Relativity and some extended theories of grav- 
itation, this method has been applied to several inter- 
esting classes of black holes coming from string theory, 
braneworlds and wormholes ^18.] . Some limitations were 
removed in Refs. fl9l l20l |. while time delay analysis was 
performed in Ref. |2ll |. 

As regards spinning black holes, the story is more com- 
plicated. Almost forty years have passed since Carter re- 
duced the geodesies equations in Kerr spacetime to first 
order equations depending on four constants of motion 
|22]. This fundamental achievement allowed a complete 
study and classification of all possible trajectories of mas- 
sive particles and photons moving around spinning black 
holes [2^. In order to visualize and study these geodesies, 
a very large amount of numerical methods has been de- 
veloped through years. In the context of gravitational 
lensing, these methods have been used to describe the 
light curve of a star orbiting a black hole 1231 and the 
apparent shape of the accretion disk |^ l25j| . Ranch 
& Blandford have proved the formation of extended 4- 
cusped caustics numerically ^26] . 

Extending the SDL methodology to axially symmetric 
black holes is not immediate and the simplicity of the 
approach may be easily lost. In Ref. [27] the SDL formula 
was recovered for light rays lying close to the equatorial 
plane of a Kerr black hole, but the coefficients of the 
formula had to be calculated numerically as functions of 
the black hole spin. Vazquez & Esteban solved the lens 
equation far from the equatorial plane for some particular 
cases IQ, but a complete analytical treatment of Kerr 
lensing is still missing. 

In this paper we make a considerable step toward this 
objective, focusing on observers lying on the equatorial 
plane and solving the general lens equation for small val- 
ues of the black hole spin. Perturbative methods allow us 
to use the Schwarzschild SDL formula as starting point 
to describe the deflection of light rays looping around 
a Kerr black hole in a completely analytical way. Our 
treatment leads to an amazingly simple description of 
all relativistic caustics and to the immediate inversion of 
lens mapping for sources near caustics. The limitation to 
the equatorial observer is motivated by the fact that the 
most important candidate black hole, Sgr A*, is likely 
to have a spin axis perpendicular to the galactic plane, 
where the solar system lies, in a first approximation. It 
is natural, then, to take advantage of this configuration 
and deal with considerably simplified equations. 

Our paper is structured as follows. Sect. 2 recalls the 
main results of Kerr geodesies. Sect. 3 explains how the 
SDL is introduced in Kerr gravitational lensing and il- 
lustrates the strategy we use to solve the geodesies equa- 
tions. Sect. 4 contains the derivation of the caustics 
order by order. In Sect. 5 we analyze the lens map close 
to the relativistic caustics, finding the additional images 



and their magnification. Sect. 6 considers the perspec- 
tives for observations in the light of what we have found, 
focusing on the study of relativistic images around Sgr 
A*. Sect. 7 summarizes the main results of the paper. 
Two appendices complement the calculations explained 
in Sect. 3 with more details. 



II. KERR GEODESICS 



In Boyer-Lindquist coordinates [2£ 
the Kerr metric reads 
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where a is the specific angular momentum of the black 
hole. All distances are measured in Schwarzschild radii 



(2A/G/C 



1). "d and (/) represent the polar and az- 



imuthal angles respectively, while x is the radial coordi- 
nate. The event horizon is a spherical surface of radius 
xu = (1 -f \J\ — Aa?')/2. In our notations, a runs from 
(Schwarzschild black hole) to 1/2 (extremal Kerr black 
hole) . 

Carter showed that the Kerr geodesies can be resolved 
in terms of first integrals of motion 22\. The final ex- 
pressions for li ghtl ike geodesies take the following form 
(following Ref. ^) 
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In these expressions, J and Q are two constants of mo- 
tion that, along with the initial condition 0o, completely 
identify the geodesic. The double signs in front of the 
integrals in Eq. remind that the integrals must be 
performed piecewise, between two consecutive values of 
X and 1? that annihilate the denominators R and O (in- 
version points). Then the sign of each piece between two 
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inversion points is chosen in such a way that all of them 
sum up with the same sign, giving a final positive contri- 
bution. 

In the context of gravitational lensing, we are inter- 
ested to those photons that come from an infinite dis- 
tance, approach the black hole reaching a minimum dis- 
tance xq and then escape back to infinity. This se- 
lects trajectories characterized by Q > 0. Moreover, 
since the roots of R represent inversion points in the 
radial motion, we have to impose that R has one non- 
degenerate positive root. This amounts to require that 
R{xq) = 0, R'{xo) 0. The limiting situation, when xq 
becomes a degenerate root, is obtained when the equa- 
tions R{xm) = 0, R'{xm) — are simultaneously fulfilled 
at some point Xm- Solving these equations w.r.t. J and 
Q we get 



4, (-3 + 2x„,) +a'^{l + 2x^) 

a(l - 2xm) 
xl, [2a^-a;„(x^-3/2)^] 
a^{x^ - 1/2)2 



(8) 
(9) 



Given a value of Xm , the quantities J„i and Qm represent 
the values of J and Q that characterize those trajectories 
that bring a photon down to the distance Xm in an infinite 
time. Asymptotically the photon keeps orbiting forever 
at a fixed distance Xm from the black hole. However, this 
orbit is unstable and small perturbations make the pho- 
ton fall into the black hole or deviate it back to infinity. 
In Schwarzschild black hole, the radius Xm of the unstable 
photon orbit is fixed to 3/2 in Schwarzschild units (the 
sphere of radius Xm is then called photon sphere) . In the 
case of Kerr black holes, the radius of the orbit depends 
on the initial orientation of the photon trajectory. In 
practice, Xm niay vary between two limiting values Xm+, 
Xm-, which respectively represent the radius of the or- 
bit described by a photon co-rotating with the black hole 
and the radius of the orbit of a counter-rotating photon 
in the equatorial plane. All intermediate values corre- 
spond to photons whose orbits are not equatorial and do 
not lie on a single plane. In order to find these limiting 



values, we have to impose — 0. So, 



and 



are found as the two largest roots of this equation. This 
is a third degree equation whose solution can be found 
exactly. However, since the successive calculations would 
soon become too cumbersome, we will already consider 
their expansions in powers of a. To describe second order 
effects in the lens equation, it is necessary to take terms 
up to the third order: 
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We see that in the limit a ^ 0, the two limiting values 
converge to the Schwarzschild photon sphere Xm = 3/2. 
When a is different from zero, Xm+ and Xm- are dis- 
tinct. The specific value of Xm in the interval [xm+, Xm-] 
uniquely fixes the amplitude of the oscillations on the 




FIG. 1: The limiting values Jm and Qm for the constants 
of motion J and Q corresponding to trajectories reaching the 
unstable orbit around the black hole asymptotically. The solid 
line is for a = Q, the dashed line is for a = 0.1 and the dotted 
line is for a — 0.2. 



equatorial plane performed by the photon along its or- 
bit. In consideration of this fact we introduce a more 
convenient parametrization, replacing a by a^: 
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(11) 



Varying the parameter ^ in the range [—1,1] we obtain 
all possible values of x,„ in the range [xm+j a;,„_], cor- 
responding to orbits with different amplitude of the os- 
cillations on the equatorial plane (a different parameter- 
ization with similar properties was used in Ref. ^Q). 
We shall see that all quantities assume very simple ex- 
pressions in terms of this parameter ^. Now, using this 
parametrization in Eqs. we can expand J„i and 

Qm to second order in a and read them as functions of 



JmiO 



'^e-2a-^!ii^ + 0(a3), (12) 
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27 



a/3 



Qrnio = ^(1 - e) + 3a^e{i - e) + ©(a^). (13) 

Notice that the presence of a in the denominators of Eqs. 
(IHl)-® allows ^ to appear at zero order already. That is 
why we needed a third order expansion for Xm- So, even 
in the Schwarzschild limit, ^ can be used to parametrize 
the photon sphere in the (J, Q) plane. 

In Fig. n^e plot the locus of points (JmTQm) when 
we vary ^ in the range [—1,1], for different values of a. 
We recall that purely prograde photons travelling on the 
equatorial plane are characterized by Q = and posi- 
tive J, while retrograde photons have negative J. Pho- 
tons with J — and Q > run on polar trajectories. 
Any geodesies characterized by J and Q outside this lo- 
cus (with Q > 0) , correspond to acceptable gravitational 
lensing trajectories. All photons with J and Q inside this 
locus are destined to fall inside the black hole. 
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There is an immediate connection between these con- 
stants and the point in the sky where the observer detects 
the photon. Throughout the paper, we consider a static 
observer at a radial Boyer-Lindquist coordinate Dq l , ly- 
ing on the equatorial plane. This restriction will keep 
all equations simple enough to be solved, while ensuring 
an exhaustive description of the expected phenomenol- 
ogy for Sgr A*. This definition has no ambiguity from a 
mathematical point of view, but needs to be linked to the 
astrophysical notion of distance from the solar system to 
Sgr A*. The current measurements of the distance to 
the Galactic center are typically based on the dynamical 
investigation of the stars orbiting around Sgr A*. The 
orbital fits are done in the context of Newtonian gravity. 
As a consequence, the current estimate of the distance 
to the Galactic center, which amounts to about 8 kpc 
0, assumes a flat background geometry. This flat 
distance makes sense as long as all scales are much larger 
than the Schwarzschild radius of the central black hole, 
which is 



-Rsch = 



2GM 



= 1.1x10' 



(14) 



for M = 3.61 X 10*5 Mq 0. Now it is evident that in 
the limit of large distances, one can simply translate any 
flat distance from the black hole, as calculated by New- 
tonian physics, into Boyer-Lindquist coordinates in the 
Kerr geometry. In fact, far from the black hole, in the 
asymptotic region, the Boyer-Lindquist coordinate coin- 
cide with the euclidean spherical coordinates centered on 
the black hole. The ambiguity in this identiflcation is of 
the order of Rsch/x {x being the distance from the black 
hole) and becomes relevant only close to the event hori- 
zon, where Newtonian physics loses any meaning. So, we 
can safely assume Dql = 8 kpc, when speaking about 
Sgr A* in any calculations. 

Then, considering only observers in the asymptotic 
region {Dql ^ 1), where the geometry is close to be 
Minkowskian, it is possible to deflne angular coordinates 
{91,62) in the observer sky. We will put the black hole 
in (0,0), and let 9i run parallel to the equatorial plane 
of the black hole while 62 will run on the perpendicular 
axis (see Fig. As Dql 3> 1, ^1 and 62 will always be 
assumed to be very small. As shown in Ref. [2^ , photons 
reaching the observer from the generic point {9i, O2) are 
characterized by the constants 



J = -e^DoL 



(15) 
(16) 



w.r.t. to the black hole. We have taken the spin axis of 
the black hole parallel to the 62 direction and we have 
considered a photon moving toward the observer. Then, 
it is immediate to pass from (J, Q) to the corresponding 
coordinates in the observer sky {61,62) and viceversa, 
apart from an ambiguity of sign in 62- 

We can use these formulae to translate the locus 
{Jm, Qm) into anew one {6i,m,62,m) in the plane (611, 02)- 
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FIG. 2: The shadow shape in the observer sky. The solid 
line is for o = 0, the dashed line is for a = 0.1 and the dotted 
line is for a = 0.2. 



This is given by 



Dol6i,„,{0 = - 
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V3 



+ 0(«')(17) 
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+ 0(a3). (18) 



and is called the shadow of the black hole. 

Fig. 13 shows the shape of the shadow in the observer 
sky for different values of a. From what we have said be- 
fore, all photons deflected by the black hole must reach 
the observer from directions {61,62) outside the shadow. 
Photons reaching the observer from the inside of the 
shadow cannot come from gravitational deflection but 
must have been generated by sources in front of the black 
hole. So, if we had enough resolution to fully image a 
black hole, we would see a black shadow with the shape 
described by Eqs. l(T7|l and ((TH|l . bordered by a luminous 
ring due to gravitational lensing of all sources around the 
black hole "q. 

In Fig. El we see that the Schwarzschild shadow is 
circular. Increasing the black hole spin a, the shadow is 
slightly distorted and gets displaced to the right, meaning 
that prograde photons (coming from the left side as seen 
from the observer) are allowed to get closer to the black 
hole, while retrograde photons (coming from the right 
side) must keep farther. 

Here, for later convenience, we are introducing and 
making use of expressions expanded to the second order 
in a. However, the exact expression for the shadow can 
be easily derived, combining Eqs. ©-(PI with Eqs. H15|) - 
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FIG. 3: The relative variation in the radial distance of the 
black hole shadow shape between the exact solution and the 
second order approximation as a function of the azimuthal 
variable ^. The solid line is for a — 0.1, the dashed line is for 
a = 0.2. 



(|16|l . Comparing the exact shadow to its second order 
approximation, we find that the latter works surprisingly 
well up to very high values of the black hole spin. In 
Fig. 131 we plot the relative error in the radial angular 
distance of the apparent shadow 9m in the approximate 
solution w.r.t the exact one as a function of the variable 
^, which follows the azimuthal angle in the observer's sky. 
Up to a < 0.28, relative variations are well under 1%, at 
a = 0.4 we find an error of 2%, while in the extremal 
case a = 0.5 the error only reaches 5%. Such errors 
must be compared to the displacement of the relativistic 
rings from the shadow, see Sect. IV, and turn out to be 
negligible for the higher order critical curves up to large 
values of the spin. 



III. KERR LENSING IN THE STRONG 
DEFLECTION LIMIT 

It is useful to introduce the following parametrization: 



9i(e,0 = ^i,™(0(l + e) 
92(e,0 = ^2,™(0(l + e) 



(19) 



Varying ^ in the range [—1, 1] and e in the range [—1, oo], 
we can obviously cover the whole upper half of the ob- 
server sky, since ^ establishes the anomaly of the light 
ray w.r.t. a reference axis in the sky (through Eqs. H17|) - 
(|18|l ) and e fixes the angular distance from the center of 
the black hole. 

In this paper, we are interested into light rays experi- 
encing very large deflections by a Kerr black hole. These 
light rays reach the observer from directions {9i, 62) very 
close to the shadow. In the parametrization H19I) . they 
are thus described by light rays with very small posi- 
tive e, while keeping ^ in the whole range [—1,1]. The 
SDL amounts to performing the integrals in the geodesies 
equations Q-©, to the lowest orders in e. 



Now, we can easily derive the values of J and Q for 
these strongly deflected photons using Eqs. (|15|1 - H16|1 : 



^(e,e)= Jrn(0(l + e) 
Q(?,e)= Q™(0(l + 2e). 



(20) 
(21) 



Substituting these expressions in Eq. Q and solving 
the equation i? = for xq, we get the closest approach 
distance as 



1 



3^/3 27 ^ ^ 



3) 



(22) 
(23) 



In general, we see that the relation between 5 and 
e depends on ^, contrarily to what happens in the 
Schwarzschild case, which, by the way, is correctly re- 
covered when a is set to zero (compare with Ref. .ITJ). 
In the resolution of the geodesies equation we will mostly 
use 5 rather than e. However, they can be immediately 
interchanged by Eq. 1231) and its inverse. 



A. Resolution strategy 

Let us introduce our gravitational lensing conflgura- 
tion. As said before, we restrict to observers on the 
equatorial plane of the black hole at radial coordinate 
Dql- We choose the zero of the azimuthal Boyer- 
Lindquist coordinate (j) on the observer position. The 
source is assumed to be static at Boyer-Lindquist coor- 
dinates {Dlsi'^s,4's)- To make contact with previous 
works, we call j — (f)s ~ tt. 

Our lens equations are provided by Eqs. Q)-® , where 
we identify (f)f — 0, (f)o = (f>s- In these equations there are 
four different integrals to solve: 



h=2 



dx 



, x^ + — aJ 
2 / z= dx 



Ji = ± / — di? 

/e 



J2 = ± 



CSC -d 



dd. 



(24) 

(25) 

(26) 
(27) 



In the radial integrals Ii and I2 we have taken the higher 
extrema to be infinite, because we assume DqltDls 3> 
1. As the two integrands go to zero as x~'^ for x —> 00, the 
relative errors committed in this approximation are of or- 
der Dq]^ and D^g respectively. Moreover, since the only 
inversion point in the radial motion is xq, the infalling 
pieces and the outgoing pieces of the integral are equal 
and we can solve the sign ambiguity considering only the 
outgoing pieces and putting a factor 2 in front of the in- 
tegrals. The radial integrals Ii and I2 can then be solved 
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using the SDL technique explained in Ref. |3- In prac- 
tice, considering photons with minimum distance very 
close to some Xm, and introducing the parametrizations 
(|20|l . (|21|l . I|22|l . we can expand the integrals in terms of 
the parameter 6, introduced before. The leading terms 
diverge logarithmically as S goes to zero, while the next- 
to-leading order terms are constants in 6. The details of 
this procedure are reported in Appendix A. Here we just 
rewrite the final results 



Ii = -fli log (5 + bi 
1-2 = -^2 log (5 + 62, 



(28) 
(29) 



where the coefficients and bi are functions of a and 
^. Their expansions to second order in a are reported in 
Appendix A. 

As regards the angular integrals (PS|) - (|?7jl . it is conve- 
nient to introduce the new variable fi = cost?. The final 
results, expanded to second order in a are functions of ^ 
and the source position /is = cosds- They are reported 
with a full derivation in Appendix B. We will just recall 
them in the following sections when we need them. 

Once all integrals are calculated, we have to solve Eqs. 
Q-© in terms of the source coordinates (7, /is), so that 
they are in the lens map form 



(30) 



Since all transformations from (7,t9s) to (7, /is) and 
from {01,02} to ((5,^) are non-singular (except for the 
points ^ = ±l)j the singularities of the Jacobian of the 
map represent gravitational lensing critical points. 

In the following sections, we will calculate the critical 
curves and the caustics of the Kerr gravitational lens or- 
der by order. Then we will describe the lens mapping in 
a neighborhood of a generic caustic, deriving the images 
position and magnification. 



IV. DERIVATION OF THE RELATIVISTIC 
CAUSTICS 

A. Zero order caustics 

Sending a to zero, we must recover the Schwarzschild 
results, i.e. that critical curves are concentric rings corre- 
sponding to point-like caustics aligned on the optical axis, 
alternatively located behind and in front of the black 
hole. Of course, as a ^ 0, all expressions are consider- 
ably simplified, and it is possible to follow calculations 
without too much effort. 

Reading all the zero-order results for the integrals in 
Appendices A and B, Eq. becomes 



-2 log (5 +21og[12(2- %/3)] = 
Ms 



= ± arcsin - 



(31) 



Defining the new variable 

^ = -21og(5 + 21og[12(2- n/3)], 
Eq. can be easily solved as 



fis = ivT-C^sinV'. 
The second lens equation ((Sj now reads 



7 — — (m — 1)71 ^ arctan ■ 



v/1 - - e 

Using Eq. to eliminate fis, we find 

7 = — (m — l)7r — arctan [^tani/;] . 



(32) 



(33) 



(34) 



(35) 



The number m appearing in this equation is the inte- 
ger part of {tp/n + 1/2) and must be interpreted as the 
number of inversions in the polar motion of the photon. 

Eqs. and represent the lens equations for a 

Schwarzschild black hole without the classical identifica- 
tion of the equatorial plane with the source-lens-observer 
plane. We can recover the results of Ref. .17] imposing 
that the motion takes place on the equatorial plane, i.e. 
setting ^ = 1. Then we have /is = (the source must co- 
herently lie on the equatorial plane) and 7 — —tp+Tr. The 
quantity ip — ir represents the deflection angle of a photon 
approaching the Schwarzschild black hole at a distance 
2^0 — Xm{^ + S). Eqs. and can also be obtained 
from the traditional planar treatment by a trivial rota- 
tion by an angle arccos^ of the reference plane around 
the optical axis using some spherical trigonometry. 

Now we can easily calculate the Jacobian of our lens 
map. We find 



and 



D 



9-0 

dj 
dip 



dS, dtp 



: sintp 



±\/l — cos V' 

tanip 
~ 1 ^2 tan2 ^ 

^ sec^ -0 
^l + S,^ tan^ %p 



dip dS, 



±- 



sin ip 

Vi^ 



(36) 
(37) 
(38) 
(39) 

(40) 



The critical curves are the solutions of the equation 
D — 0, which, in our case, simply gives 



ipcr — kn. 



(41) 



The critical ip does not depend on ^. Recalling that ip 
is just a function of S expressed by Eq. we have 



12(2- V3)e 



-kn/2 



(42) 
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Switching to e by Eq. 123|) . we have 



(43) 



Then, recaUing the meaning of e by Eq. (|19|l and taking 
6'i_m and 6*2,™ from Eqs. (|17|l - p8|) . we finally find 



DoLdl,cr{£.) 
DoLd2,cr{£.) 



(44) 



that is a series of rings parameterized by ^, slightly larger 
than the shadow of the black hole, in full agreement with 
Refs. H mill 113. The critical curves are labelled 
by the number k. We shall refer to k as the critical 
curve order or the caustic order, when we consider the 
corresponding caustic. 

Coming to the caustics, inserting Eq. (|41|l in the lens 
equations Eqs. H33() - (|35|l and noting that the number of 
inversions in polar motion is m = k, we find 



^J■s 



0, 7 = -(fc-l)7r. 



(45) 



The caustics are points aligned on the optical axis. For 
odd k they are behind the black hole, while for even k 
they are in front of the black hole (retro-lensing caus- 
tics). In the Schwarzschild limit, the number of loops 
performed by photons forming critical images of order k 
is just (fc — l)/2. However, this is not exactly true for 
spinning black holes, as we shall see in the next section. 
So, it is better to think of the order of the critical curve 
k as the number of inversions in the polar motion done 
by photons associated to it. 

Of course, the first caustic for fc = 1 is the classical 
weak field limit one, which is outside the range of the 
SDL approximation, so we cannot expect to recover the 
Einstein ring radius putting k = 1 into Eq. (|44|l (the 
first caustic is no longer described even by the weak field 
approximation if the source is close to the black hole). 
However, as shown in Refs. 0,01, the SDL approxima- 
tion works better and better for higher k, starting from 
the first retro-lensing caustic in fc = 2. It is these caustics 
that we are going to study in the following sections. In 
particular, we will find out how they are displaced and 
deformed by the black hole spin, obtaining a full analytic 
description of their shape. 



B. First order caustics 

Up to now we have just re-obtained all the 
Schwarzschild black hole results in a more complicated 
form, starting from Kerr geodesies equations and send- 
ing a back to zero. Now, we shall introduce first order 
corrections to our lens equations, re-deriving the critical 
curves and the caustics. We anticipate that the caus- 
tics get displaced from the optical axis in the azimuthal 
direction, though remaining point-like. 

Using the first order terms of the radial and angular 
integrals from the Appendices A and B, we can add the 



terms proportional to a in the equations H33|l . (|35|l . The 

inversion of the /is equation can be easily performed order 
by order in a. Then, repeating the same steps of the 
previous section, we find 



4a^ 



= ±^1 -"^^ sinV-i -^y/l-e sinV' 
3v 3 



(46) 



l)7r — arctan tani/;] 



4a 
37^ 



tan?/) 



1- (l~,^2)sin^^ 



tani/>-|-3%/31og(2\/3-3) 



(47) 



Note that for ?/' close to imr + 7r/2 and ^ close to zero, 
the first order correction to fis may bring it to absolute 
values larger than 1. As fig is the cosine of the polar 
angle, these values are unphysical. This inconsistency 
comes out because, when we solve for fis order by order, 
we expand the arcsin function in points very close to 
the extrema of its definition range, where the arcsin is 
not analytic. Then, the linear approximation obviously 
takes us out of the interval [—1,1]. The dangerous values 
of and ^ correspond to nearly polar trajectories where 
the final direction is very close to one of the two poles. 
However, as we shall see, the highest magnification for 
the relativistic images is obtained when the source is close 
to a caustic. Luckily, the caustics lie at ~ for an 
equatorial observer, so that we will always work very far 
from these dangerous regions. This danger will become 
effective for very high order caustics, which may become 
very large, as we shall see in the next subsection. 

Now we can calculate the derivatives of the lens equa- 
tion as before. The Jacobian reads 

sinV' SaS, \ 



D = ±- 



1 



3%/3; ' 



(48) 



which tells us that the critical curves are still described 
by Eq. (|41|l even at first order. This means that there is 
no correction to the critical 6 H42I) of the previous section. 
The fact that we do not get any corrections to 5 does not 
mean that the shape of the critical curves is not altered 
by the black hole spin at first order. Indeed going back 
from 6 to £ we get a first order correction, according 
to Eq. H23(l . Moreover, the shadow shape is modified 



according to Eqs. 
critical curves are 



||T7J|-||TH1). Then, to first order, the 



DoLSl.criO — 
DoL62,cr{0 = 



3\/3g 
2. 

+2a [1 
±2ae, 



(1 + ecr) 



1 + Cnr 



(49) 



with ecr still given by Eq. H43() . 

Also the caustics are modified. In fact, plugging ipcr = 
k-n- into Eqs. ignil-lEZI), we find 

Ms = (50) 
kir 



7 — —{k — l)n — 4a 



3\/3 



log(2V3 - 3) 



(51) 
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The caustics are still fully confined to the equatorial 
plane, they are still point-like, but they are displaced 
from the optical axis. The displacement is negative, 
which means that the caustics drift clockwise if we look 
at the Kerr black hole from the northern pole. So, we 
can say that if a source lies on a caustic point of order 
fc, prograde light rays perform more than (fc— l)/2 loops 
while retrograde light rays perform less than (fc — l)/2 
loops. The number of inversions in the polar motion is 
still k. Higher order caustics are more displaced, because 
of the k dependence in Eq. (|51|) . Of course, as said be- 
fore, this formula correctly describes all caustics except 
for the weak field one, corresponding to A: = 1. Going to 
second order in a we will describe the full shape of the 
caustics. 



as we already know the zero order critical curve and we 
know that there is no correction at the first order. Then 
we easily get rid of all the trigonometric functions and 
the final Jacobian reads 



9 Sip- (92 - 24V3 - 15A:7r)(l - S,^) 



(55) 

The (—!)'" is a consequence of the expansions of the 
trigonometric functions, while the double sign inherited 
by the Jacobian at all orders depends on the fact that 
the (e, ^) parametrization only covers half of the observer 
sky and we are forced to introduce a double sign in the 
equation for ^s- 

The equation D = gives the second order correction 
to the critical ijj in a very simple form 



C. Second Order Caustics 

At first order in a the caustics still remain point-like, 
while it is known that they get a finite extension when a is 
different from zero 26, 27]. So, it is necessary to proceed 
to second order. To the right hand side of Eqs. lP?i|l - l|T7| 
we have to add the following quantities respectively: 



a X 



1273(1-^2)^^ cosV 
-V3(l -C2)sin3V' 



(52) 



2a- 
27l 



2 (1 - ^^)(23C2 _ 4) arctan(Ct) 



t f 8 l + t^ 



a'i [(1 -f 



(53) 



where t = tan?/' and '0„ = 5?/) H- 8V3 — 20. 

Note that a term directly proportional to ip appears 
in Eq. (|52|l through ?/;„. This is another danger for the 
approximation, since for very large "0, i.e. photons per- 

forming several loops around the black hole, 5^s may 
become even larger than 1. This break-down sets the 
true limit to the applicability range of the perturbative 
expansion in a, which gets smaller and smaller for pho- 
tons making many loops. However, the brightest images 
are formed by photons associated to critical images of low 
order. For these images, as we shall see, the range of ap- 
plicability of the perturbative expansion is considerably 
large. 

Now let us find out the corrections to the critical 
curves. Once we have calculated all derivatives of the 
lens equations and written the Jacobian to the second 
order in a, we set 



tp = kn + a Sijj, 



(54) 



92_24V3-15fc7r ^ 

H= g (1-C)- (56) 

Now we can remount the complete second order ex- 
pansion of the critical curves, which reads 



DoLei,cr[0 



-2a [l + e„{l-e)\ 
-igU^-2+[ffc7r(l-a 
^116 _ 12^3 +(ip_i2V3) e] 

3V3 



a-{e+[^kn{l~e) 



± 



V3_ 



+ 12VS - 54 + (ip - 12V3) e] ecr} 

(57) 

Here, again, we have kir terms which become large for 
higher order critical curves. 

Finally, let us calculate the caustics at the second order 
in a. Plugging Eq. with ifS^ into the lens equations, 
we find 



±(-l)V,(l-e^)3/2 
-(fc-l)7r-4a ^'"^ 



3\/3 



log(2x/3 - 3) 



-ret 



where we define 



7'c = —a 
9 



'(5fc7r-|-8\/3-36). 



(58) 



(59) 



(60) 



The analytical expressions of the Kerr black hole caus- 
tics, given by Eqs. (|58|I - H59() to the second order in the 
black hole spin a, represent a major achievement of this 
paper. Before discussing their shape and all the physical 
implications, it is a good idea to test our formulae by 
comparing them to the results obtained by alternative 
methods. In Ref. [23|, the intersections of the caustics 
with the equatorial plane were found using the SDL ap- 
proximation only, without any limitation for the black 
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a 



-0.4 -0.2 0.2 0.4 



FIG. 4: Comparison between the intersections of the caustics 
with the equatorial plane as calculated in Ref. (27| without 
the perturbative approximation for the black hole spin (solid 
lines) and the ones calculated in the present paper (dashed 
lines). The plot refers to the caustics of order 2 < fc < 7. 



hole spin. The first test is to analytically expand the for- 
mulae of Ref. Il^l to the second order in a, without us- 
ing any numerical integrations. Indeed, we get the same 
expressions as in Eq. H59|l . evaluated for ^ = ±1. Fur- 
thermore, we can draw in the same plot the intersections 
of the caustics with the equatorial plane as calculated 
in this paper along with those calculated in Ref. [2^ . 
Rather than making two separate plots for prograde and 
retrograde photons, we can make a unique plot, letting 
a vary in the range [—0.5, 0.5] and keeping the values of 
Eq. ^ for C = 1- So, the left side of the plot (a < 0) 
represents the intersections for retrograde photons and 
the right side (a > 0) represents the intersections of pro- 
grade photons. We see in Fig. ^that the second order 
approximations (dashed lines) follow the exact expres- 
sions of Ref. very accurately. We can estimate that 
for lower order caustics the perturbative approximation 
works up to a ~ 0.3, while for the last caustic in the 
plot {k = 7) we have to stop at a ~ 0.1. In any case, 
the validity range is impressively large, reaching values 
of the black hole spin comparable to the extremal case. 
This encourages us to make extensive and confident use 
of the second order approximation for a full description 
of Kerr lensing phenomenology. 

Now, let us discuss the shape and the extension of Kerr 
lensing caustics. Looking at Eq. (|58|l . it is interesting to 
note that the upper half of the critical curve is mapped in 
the lower half of the caustic for odd k, while it is mapped 
in the upper half if k is even. As already found numer- 
ically in Ref. |2^, the caustics have the characteristic 
astroid shape shown in Fig. [S] common to all tangential 
caustics after the breaking of the axial symmetry. The 
four cusps are in ^ = ±1 and ^ = choosing different 
signs for ^j. 

The caustics have the same extension in 7 and /i. We 
recall that 7 is just the azimuthal angle of the Boyer- 
Lindquist coordinates taken from the reference axis start- 



:i 




FIG. 5: The typical caustic in Kerr gravitational lensing has 
the astroid shape and the same angular extension Tc (given 
by Eq. 1601 ') along the azimuthal and the polar direction. 





a=0.01 


a=0.05 


a=0.1 


a=0.2 


A72 


0.018 


0.088 


0.18 


0.35 


A73 


0.042 


0.21 


0.42 


0.84 


A74 


0.066 


0.33 


0.66 


1.3 


A75 


0.09 


0.45 


0.9 


1.8 


A76 


0.11 


0.57 


1.14 


2.3 


A77 


0.14 


0.69 


1.39 


2.8 



TABLE I: Drift (in radians) of the caustics of order k with 
2 < fc < 7, for different values of a. 



ing from the black hole and going in the direction oppo- 
site to the observer, while /i is the cosine of the polar 
angle d. As said before, we can trust our results as long 
as the perturbative terms remain small. In this regime, 
/i ~ § — ^- Then, we deduce that the caustics have the 
same extension in the azimuthal and in the polar direc- 
tion. 

The extension is parameterized by the quantity Tc, 
which is the semi- axis of the caustic. We see that it 
grows with the black hole spin a and with the caustic or- 
der k. Thus, for higher caustic orders, the perturbative 
approximation fails for smaller and smaller values of a, 
while it remains good for lower order caustics. This was 
already noted while commenting Fig. ^ 

The drift from the optical axis of the caustic is given 
by the first order term in Eq. (|59|l . Indicating it by A7, 
in Tab. 1 we report the obtained values for the first 6 
relativistic caustics, starting from fc = 2, for values of 
the black hole spin going up to a — 0.2. We see that 
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a=0.01 


a=0.05 


a=0.1 


a=0.2 




2 


0.00021 


0.0052 


0.021 


0.082 


Tc 


3 


0.00056 


0.014 


0.056 


0.22 


Tc 


4 


0.0009 


0.023 


0.09 


0.36 


rc 


5 


0.0013 


0.031 


0.13 


0.5 


Tc 


6 


0.0016 


0.04 


0.16 


0.64 


Vc 


7 


0.002 


0.05 


0.2 


0.78 





TABLE II: Radius (in radians) of the caustics of order k with 
2 < k < 7, for different values of a. 



the drift may become very large, amounting to several 
radians already for the 5th order caustic, while still in 
a regime where the perturbative solution is valid. For 
higher order caustics, the number of loops may be signif- 
icantly different from the planar orbit result (fc — l)/2. 
Another important consideration is that we do not need 
perfect alignment between source, lens and observer to 
have gravitational lensing. To enhance the images asso- 
ciated with a critical curve of order k, the source must 
align with the corresponding order k caustic, which may 
be well off the optical axis. Moreover, the relativistic 
images are not enhanced all at at the same time, since 
caustics of different order move far away each other. 

In Table 2 we report the radii of the first 6 relativistic 
caustics for different values of a. The extension of the 
caustics is of second order in a and thus remains much 
smaller than the drift, reaching some tenth of radians in 
the perturbative regime. Outside of this regime, it is difh- 
cult to know what would happen to higher order caustics. 
Would their vertical extension saturate before reaching 
the poles or would they wrap around the pole? Would 
they meet each other and make transitions to more com- 
plicated structures? The answers to these questions need 
further research, both analytically and numerically. We 
just want to remark that the finite extension of relativis- 
tic caustics is of striking importance for phenomenology, 
as will be clear in the next sections. 



V. GRAVITATIONAL LENSING NEAR 
CAUSTICS 

The description of the caustics is the fundamental step 
for a full description of gravitational lensing. In this sec- 
tion we will give a complete analytic resolution of the 
Kerr lens equation for sources close to relativistic caus- 
tics. 

The starting point is the second order lens equations, 
built adding to (gHl and to igTI). Let us consider 
a source whose distance from the fc-th order caustic is 
of order (thus being comparable to the caustic size). 
Then its position can be expressed in the following way 



Ms 
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icau 



(61) 
(62) 



where ^caL is the caustic position at the first order in a, 
expressed by Eq. ifCT) . Correspondingly, the images as- 
sociated to the critical curve of order k will be enhanced. 
They will be described by 



-0 = fcTT + a^Stp. 



(63) 



Substituting Eqs. (|61|I - H63I) in the lens equations, the 
zero and first order terms cancel out and we are only left 
with the second order terms 



^7 ■ 



(l-C^)(5fc7r + 8\/3-20) 



32 



+ {3 - f ){5kn + SVS - 20) 



(64) 



(65) 



This is the Kerr lens equation close to the caustic of 
order k. S is just a sign which takes into account the 
fact that the (V'jC) parametrization only covers half of 
the observer sky. So, S = +1 for the upper half of the 
observer sky and S = —1 for the lower half. We can 



easily check that the Jacobian 
Jacobian of this lens equation 



is just a times the 



D 



(66) 



The surprisingly simple form of the lens equation en- 
courages its analytical resolution. The 5^ equation can 
be easily solved for dif): 



A / '^^ ^ 1 



32 - (3 - $2) (5fc7r -I- SVS - 20)1 . (67) 



Plugging this expression into the (5/is equation, we can 
write it in the form 



(68) 



where Xc — Tc/o? and rc is the semi-axis of the caustic as 
defined by Eq. I|60() . Squaring both sides we get a fourth 
order equation for ^ 

xli'^ + 2xJ-r^^ + (S-f^ + S^il ~ xDi' - 2xc5ii - S-f^ = 0. 

(69) 

The real solutions of this equation are images for a 
source in {-^cau + , a^Sfis)- It is easy to check that we 
have two images if the source is outside the caustic and 
four images if it is inside. 

Once we have found the solutions of the squared equa- 
tion, we have to go back to the original equation H68|l . 
Each root of Eq. (|69|l satisfies Eq. (|68|l only with one 
choice of S. This determines the half-sky where the im- 
age appears. It is the upper half if 5 = 1 and the lower 
half if S* = — 1. Then, we can easily calculate the value 
of Sip for each image through (|57|l and then go back to e 
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by Eqs. ^ and Finally we can write the images 

as 

DolOi = + ecr) + 2a[l + e,,(l - C')] 



[12-6^2 + ecr(76- 82^2_|_27^^)] (70) 



DoL02 =^-^Vl^ 



(1 + £cr) + ^"^^'^ 



(71) 



with ^ and (^V" solving Eqs. 

In the particular case Sfx^ — (source on the equa- 
torial plane), the solutions are ^ = ±1 and ^ — Sj/xc 
(double root). The first two solutions are two images 
staying on the equatorial plane, which are physical for 
any value of 6'^. The other two are acceptable only if 
|<57| < Xc because ^ is defined in the range [—1,1]. This is 
in agreement with the fact that Xc represents the caustic 
semi-axis. These two images form symmetrically w.r.t. 
the equatorial plane, grazing the critical curve from the 
outside. When these additional images are present, the 
former two are inside the critical curve, while when they 
are absent the remaining images are one inside and the 
other outside of the critical curve so that global parity is 
conserved in caustic crossing. 

On the other hand, if Sj — 0, we have the solutions 
^ = (double root) and ^ = ± ^/T~{5ii7/xcf . The first 
two form very close to the polar direction on opposite 
sides of the critical curve, while the last two arc real only 
for \Sfj,s\ < Xc- They form symmetrically w.r.t. the polar 
direction and graze the critical curve from the inside. As 
before, things work in such a way that global parity is 
conserved. 

Now let us make an example with a physical source and 
a physical black hole. Sgr A* has a mass of 3.61 x 10^ so- 
lar masses. Let us suppose that its spin is a = 0.02 (Liu 
& Melia estimate a — 0.044 [s^l, but different methods 
point to different values). Then we are able to calculate 
the caustic positions and shapes. As a source, consider 
a star with a radius Rs = Bi?© at 200 AU from Sgr A*. 
This is the order of magnitude of the periapse distance 
of the observed stars orbiting Sgr A*, like S2 or S14 0. 
In Fig. we show what we would see if this star ap- 
proaches the first retro-lensing caustic. The position of 
this caustic is in 7 = 177.98°, so slightly displaced from 
the optical axis. This means that the source should be 
almost in front of the observer, very close to the optical 
axis. On the left panels of Fig. we show several po- 
sitions of the source relative to the caustic, as seen by 
the black hole. Notice that with these values of a, source 
radius and distance the size of the source is comparable 
to that of the caustic. On the right panels, we show the 
corresponding images and the shadow in dashed style. 
We see that when the source is far from the caustic (top 
panels) there are two small images. The bigger one is 
below the black hole if the source is above the equatorial 
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FIG. 6; Formation of the images for a source approaching 
the first restrolensing caustic (fc = 2). On the left we show 
several positions for a source and on the right we have the 
corresponding images around the shadow (in dashed style). 
The thickness of the images has been exaggerated to make 
them more evident. 



plane (we recall that this is normal in a retro-lensing sit- 
uation). When the source approaches the caustic (second 
row panels), the two images do not lie on opposite sides 
w.r.t the black hole. This distortion is a consequence of 
the axial symmetry breaking. When the source enters 
the caustic, two more images form (third and fourth row 
panels). If the source orbits with a velocity of the order 
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FIG. 7; Formation of the images for a source approaching 
the first relativistic standard lensing caustic (fc = 3). On the 
left we show several positions for a source and on the right we 
have the corresponding images around the shadow (in dashed 
style). The thickness of the images has been exaggerated to 
make them more evident. 



of the circular orbit speed -y/ GM/Dls, the whole caus- 
tic crossing takes 3.4 hours, much longer than the typ- 
ical times of the primary caustic crossing, which takes 
just few seconds [2g. Furthermore, since the higher or- 
der caustics are much more extended, the probability of 
caustic crossing is much higher. 

In Fig. [7| we have shown the case where the same 
source approaches the first relativistic standard lensing 
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FIG. 8: A zoom very close to the shadow border (in dashed 
style), showing at the same time the images of two sources, 
one being in the k — 2 caustic (outer tangential arc) and the 
other being in the fc = 3 caustic (inner tangential arc). 



caustic (fc ~ 3), which now is displaced to 7 = 4.8° on 
the right of the black hole. As this caustic is larger, 
the source now looks smaller compared to the caustic, 
as we see in left panels. When the source is far from the 
caustic (top panels), there are two images, the bigger one 
being on the same side of the source (standard lensing 
situation). As the sources approaches, the images and 
the black hole arc no longer on the same line (second row 
panels), then formation of two new images takes place 
(third and fourth row panels). In this case, the caustic 
crossing takes 9.2 hours, for a source velocity equal to 
the circular orbit velocity. 

Fig. IHl zooms on two images generated by different 
sources in the caustics k — 2 and fc = 3. The first source 
generates the outer tangential arc while the second source 
generates the inner tangential arc. This is because the 
more loops the photons perform, the closer they get to 
the black hole. Then higher order images appear closer 
and closer to the shadow. What astonishes of this picture 
is the tremendously small thickness of the arcs. These 
may be greatly elongated and even form full rings if the 
source is larger than the caustic, but their radial size is 
really very small. In Figs. 16171 we had to exaggerate the 
thickness in order to show them in a more evident way. 
The next section will be devoted to the calculation of the 
length and the thickness of the arcs, i.e. the magnifica- 
tion of the images. 
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A. Magnification 

In standard weak field gravitational lensing the magni- 
fication is the ratio between the angular area of an image 
and the angular area of the source if no lensing occurred. 
This does not necessarily make sense when high deflec- 
tion takes place, since the side of the source seen by the 
black hole is generally different from the side seen by the 
observer. Then, if the source does not emit isotropically, 
a magnification calculated in the standard way would not 
give the correct ratio between the brightness of the image 
and that of the source. For example, in the retrolensing 
situation, the source is in front of the black hole. So 
the photons going toward the black hole leave the source 
from the side opposite to the one seen by the observer. 
For simplicity, in this section we shall assume that the 
source emits isotropically. The formulae can be easily 
corrected in the case this does not happen. The isotropic 
emission hypothesis ensures that the source as seen by 
the observer is simply D'j^s/^OS smaller than as seen by 
an observer in the black hole position in the absence of 
the lens. 

The angular area of the image in the observer 
sky is simply ddid92- The angular area of the 
source in the black hole Boyer-Lindquist coordinates is 
I smi!}sd'yd'ds \ = djdfis, when the source is very far from 
the black hole. Then, the magnification matrix is just 
the Jacobian matrix of the lens map in the form 



17-7(^1,^2) 
I = Ms (6*1, ^2) 

Yet, we have Eqs. H64(l - ()f)5(l in the form 
and Eqs. in the form 

\ 02 = ^2^,0 



(72) 



(73) 



(74) 



the signs. We will not write them here, but we shall give 
the two eigenvalues of the Jacobian matrix 



Xr 



2DoL 
2DoL 



3V3(1 



+ Ecr) 



-Dn 



(77) 
(78) 



with Dq being proportional to the Jacobian studied in 
the former section (|55|l 



Do = ±^l-eD = 



(-l)fe Q' 

9 



9(5?^-(92 - 24^/3-15/^7^)(l-^2)l ^ (^^9) 



Since Ccr is fixed by the caustic order fc, is always 
positive, while At vanishes whenever Dq does. This con- 
dition is fulfilled when Eq. (|56|l holds, i.e. on critical 
images. It is possible to show that the two eigenvectors 
associated to A^ and At respectively become radial and 
tangential in the limit a — > 0. So, when the source is 
close to a caustic, all images are elongated in a direction 
nearly tangential to the critical curve, as already noticed 
in the previous subsection. We shall call A^ and At radial 
and tangential eigenvalues respectively, though they are 
such only in the limit a — + 0, actually. 

Finally, we can write the radial and tangential mag- 
nification of the images w.r.t. the source as seen by the 
observer. These are 



Dls \ 



Dos 



-cr 

1 _ Dos 3%/3(l + tcr) 



Dls I At I Dls 2Dol\D, 



(80) 
(81) 



Of course, the total magnification is /i = fJ-rfJ-t- A 
good check is to reduce this formula in the Schwarzschild 
limit to compare with Refs. |l4|, l2(l and in the 
equatorial limit, to compare with Ref. [23. The first 
limit is obtained sending a to zero keeping the source 
position {a'^Sj, a?5^,s) fixed. This is equivalent to put Xc 
to zero in Eqs. Ht)9|) . Ht)7|) . Then the images are in 



Then we can find the Jacobian matrix of Eq. H72(l as 



9(01 



where we have used the matrix notation 

dyi dyi 



dy2 dy2 
dxi dx2 



(75) 



(76) 



and we have noted that dj — a?d{5^) and d/is = 
a?d{6pLs). 

Then writing the explicit expression of the elements of 
J is straightforward, once we correctly take care of all 



5i}= - 



mill 



T V^7^ + ^- 



9(^7' + ,5^2) 
Substituting in Eq. (|81|) . we find 

27ecr(l + ecr) 



l^Sch 



^os 



(82) 
(83) 



(84) 



Identifying ^y&y^ + S]l^ with the misalignment of the 
source with the point-like caustic position, we exactly 
find the magnification of Refs. |12.l20l| . 

The equatorial limit is recovered when 6 fig = 0. Then 
we have two equatorial images plus two non-equatorial 
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images if the source is inside the caustic. The two equa- 
torial images are described by 



e= ±1 
Sip = — 7'c T Sj. 

Inserting these values in Eq. (|81|l . we find 



Meg 



^os 



27ecr(l + e-cr) 



(85) 
(86) 



(87) 



which is the leading term close to the equatorial cusps as 
found in Ref. ^ . 

The magnification of relativistic images is usually very 
low. This is expected by the fact that a very small pertur- 
bation in the photon trajectory may completely change 
its final direction. Referring to a source at lOOAU from 
Sgr A*, by Eq. (|80|l we find a radial magnification 



^1^^^ = 5.3 X 10 
nl?^ ^ 2.3 X 10"^ 



/ lOOA[/ \ 
V Dls ) 



(88) 
(89) 



for the relativistic images of order 2 and 3 respectively. 
These very low values, which are independent (at low- 
est order) of the distance of the source from the caustic, 
justify the very thin arcs of Fig. |S1 

The tangential magnification diverges when the source 
crosses a caustic. However, at most the images may 
merge to form a full Einstein ring, which gives the maxi- 
mal tangential magnification. Dividing the angular area 
of a circular corona of radius 3-\/3(l -I- ecr)/{2DoL) and 
thickness iir{'2Rs)/Dos by the angular area of the source 
'k{Rs / DosY' , we obtain the maximal total magnification 



(2) 
f^max 



f^max 



1 Dls 27e,.(l + 

^cr J 



Rs Dls 
8.7 X 10"^ 



L»2 



Rs_ 

Rq 



( Dls 
\imAU 



1 Dls27e^r{^ + e„) 
Rs Dls Dq^ 

3.8xl0-f^y' (^^ 



'LS 



(90) 



(91) 



In Fig. I^lwe show a map of the tangential magnifica- 
tion centered on the first retro-lensing caustic (the caustic 
order is fc = 2) for a = 0.02. To get the total magnifica- 
tion map, it is necessary to multiply the values on this 
map by the radial magnification 1)88(1 , which is practically 
constant in a the neighborhood of the caustic. It is inter- 
esting to note that the tangential magnification is greater 
than one in a large region surrounding the caustic. This 
means that the cross section for the generation of large 
relativistic tangential arcs is quite high. This is a very 
important characteristic of relativistic images that we are 
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FIG. 9: Tangential magnification map centered on the k = 
2 caustic (the first retro-lensing caustic) for a = 0.02 and 
Dls = IQQAU. 



going to exploit in the next section when discussing their 
observability. 

For a source at the center of the caustic, we have an 
Einstein cross where all images have the same magnifi- 
cation. Then the total magnification takes a very simple 
expression 



_ D^S 27ecr(l + Ecr) 
^LS^OL 'c 



(92) 



where Tc is the semi-axis of the caustic as defined by Eq. 
(EOJl. 



VI. PERSPECTIVES FOR OBSERVATIONS 

The relativistic images appear just outside the shadow 
of the black hole. In order to distinguish them, we need a 
resolution of the order of /xas. The present world record 
has been achieved with Very Long Baseline Interferom- 
etry (VLBI) in the mm band and amounts to 18 /xas 
f^. However, in this band, there are no good compact 
candidate sources around Sgr A* for gravitational lens- 
ing. More interesting are the infrared and especially the 
X-ray band. 



A. Infrared band 

In the infrared K-band, centered at A = 2.2 /^m, the 
exinction by interstellar dust allows good observations of 
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the stellar environment around Sgr A* Many stars 
have been detected and followed during their orbital mo- 
tion around Sgr A*, providing the best dynamical con- 
straints on its mass distribution. Surprisingly, these stars 
are of early spectral types, leaving open the question 
on the presence of such young stars in the galactic cen- 
ter. In the K-band, these stars have magnitudes between 
rriK = 13 and ttik = 16. 

As regards the angular resolution in the K-band, 
the VLT units can be combined to perform in- 
terferometry observations with an equivalent base- 
line of 200 m and a maximal angular resolu- 
tion of 2.2 mas (http://www.eso.org/projects/vlti). 
Some space missions performing nulling interfer- 
ometry (TPF, http://www.terrestrial-planet-finder.com 
DARWIN, http://ast.star.rl.ac.uk/darwin) should be 
launched in the near future. According to the mission 
designs, some spacecraft should fly in formation at dis- 
tances of the order of tens of meters. A futuristic de- 
velopment of such idea might lead to much higher res- 
olutions. The baseline needed for l^aa resolution is of 
the order of hundreds of kilometers. High precision for- 
mation flying may be achieved by laser ranging and mi- 
crothrusters in the wake of what is being studied for LISA 
(http://lisa.nasa.gov), where the distance between the 
spacecraft is 5 million km. 

Some of the stars around Sgr A* may cross some caus- 
tics and generate bright relativistic images. However, 
they would be embedded in the flux coming from Sgr 
A* environment. In the quiescent state, Sgr A* flux in 
the K-band should have rriK > 18.8 The infrared 

emission of Sgr A* is believed to originate in the inner 
10 Schwarzschild radii of the black hole [s^, with the 
Schwarzschild radius being 



2GM 



■ch 



1.1x10^" 



(93) 



Then it is necessary to establish whether relativistic im- 
ages can overcome the background flux. As conservation 
of surface brightness holds in gravitational lensing, the 
relativistic images must have the same surface brightness 
of the original source as seen by the observer, 



Is = 



Ls 



(94) 



where Ls the intrinsic luminosity of the source, Dqs — 8 
kpc is the distance to the source and fls is the angular 
area in the sky subtended by the source, which is fig ~ 
tt{RsY / D'^g . For a source of 10 solar radii with nriK = 
15, the surface brightness of the relativistic images is four 
order of magnitudes larger than that of Sgr A*. 

However, relativistic images have a tiny angular area 
and their contribution to the number of photons collected 
by a pixel in a CCD detector may be very small. To get 
an idea of this fact, let us consider a CCD detector where 
every pixel collects energy flux from an angular area of 
size If a relativistic image, in the form of a tangential 
arc of angular thickness ujarc = fJri'^Rs)/ DoS: lies on the 



pixel area, the flux received is S' oc Is{ujpXuJarc)- On the 
other hand, the noise coming from the environment of 
Sgr A* to the single pixel is iV oc IsgrUJp. Then, the 
signal-to-noise ratio for a single pixel is 



5^ 

N 



Is 



I Sgr 



2.2 



1/ias 



Rs 

10 Rr. 



(95) 

(l^)"^-(-) 



i.e. a pixel with a tangential arc receives only twice more 
K-band photons than other pixels. The S/N could be im- 
proved by taking a smaller pixel size or stars with lower 
radii and higher brightness, but we cannot go very far. 
Taking into account absorption by the matter surround- 
ing the black hole that would surely take place and fluc- 
tuations in the surface brightness of Sgr A*, we doubt 
that relativistic images of stellar sources can be actually 
detected in the K-band by present or near future tech- 
nology. 



B. X-ray band 

In the X-rays, Chandra is leading very impor- 
tant observations discovering the physics of high en- 
ergy electromagnetic sources in the central region of 
the Galaxy [33 with a resolution of the order of 
0.5 arcseco nds. The space mission project MAXIM 
(http://maxim.gsfc.nasa.gov) will represent a major leap 
toward high resolution, reaching the striking resolution 
of 0.1 /ias. With such observational facility, a complete 
and detailed imaging of the black hole will be possible. 

It is very interesting to consider that Sgr A* luminosity 
in the 2-10 keV band is 2x10'^^ ergs s~^, which is much 
lower than expected if the black hole were accreting at 



the Eddington rate, i.e L = 3x10 ergs s~ [35j- There 
are several models for Sgr A* accretion. Models based 
on Bondi accretion (spherically symmetric inflow) predict 
the X-ray emission to be created in a region of the order 
of lO^Rsch IjS^- Models based on advection dominated 
flow [33 predict the emission to be dominated by cooler 
gas at larger radii, of the order of 10'*i?sc/i- 

Indeed, many X-ray sources have been detected in the 
neighbourhood of Sgr A*, with a luminosity comparable 
or even slightly higher than the supermassive black hole 
[36| . These sources are probably Low Mass X-ray Bina- 
ries (LMXB) which seem very numerous in the galactic 
center. The situation seems really appealing, since we 
have a population of bright compact sources, with possi- 
bly poor contamination from the intrinsic luminosity of 
Sgr A*. It is believed that the most of the X-ray emission 
from a LMXB comes from a region of tens of kilometers. 
Then we can assume Rs — 10^ km, with an X-ray flux 
of Ls ~ Lsgr = 2x10"^'^ crgs s~^. If we consider an emit- 
ting region of lO^Rsch for Sgr A*, there are 14 orders of 
magnitude between the surface brightness of an LMXB 
and the surface brightness of the X-ray emission of the 
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supermassive black hole environment. Then, the signal- 
to-noise ratio for a single pixel imaging a tangential arc, 
Eq. is 

^=o.9xio«fiy' (^yi^y 

N \l^as J \100kmJ \ 100 AU J 

(97) 

i.e. the signal in a pixel touched by a tangential arc 
is nearly 6 order of magnitudes higher than the noise 
from Sgr A* environment for a detector with the accuracy 
of 1 //as. The contamination from Sgr A* environment 
seems to be completely under control. Relativistic grav- 
itational self-lensing of Sgr A* would just produce rela- 
tivistic Einstein rings with the same surface brightness 
of Sgr A*. The only serious danger for photons coming 
from outside Sgr A* and deflected by the central black 
hole is absorption by the matter surrounding it. How- 
ever, even if some absorption certainly occurs, it seems 
difficult to fill a gap of so many orders of magnitudes 
between the surface brightness of the relativistic images 
and that of Sgr A* without affecting the luminosity of 
Sgr A* as well. 

Of course, the idea firstly proposed by Ranch & Bland- 
ford that some X-ray flares may be explained by lensing 
of nearby sources is fully plausible in this scenario p^. 



VII. CONCLUSIONS 

In this paper we have made an analytical treatment 
of gravitational lensing by Kerr black holes in the strong 
deflection limit. In order to achieve our objective we have 
made three approximations. 

The first one is the strong deflection limit approxima- 
tion for all radial integrals. This is just an expansion 
of the elliptic integrals that result from integrations over 
the full radial motion of a photon. Restricting to pho- 
tons suffering a very large deflection, for all the radial 
integrals we have only kept the leading term diverging 
as log e and the constant term (e being the separation 
between the image and the shadow border as seen by the 
observer). As shown in several articles, this limit gives a 
very good approximation starting from photons deflected 
by an angle of order tt P,0,[ll|i|]. 

The second approximation has been to consider only 
small values of the black hole spin a. This has allowed 
us to take the Schwarzschild gravitational lensing as a 
starting point for the derivation of the corrections due 
to the presence of an intrinsic angular momentum of the 
black hole. As far as we could compare our results with 
available exact ones, we have verified a considerably wide 
applicability range of our approximation. For the first 
relativistic images, we can safely apply our treatment up 
to a = 0.1 {a — 0.5 being the extremal Kerr black hole 
in our units). 

The third restriction concerns the position of source 
and observer. Besides the limitation to far sources and 
observers {Dls, Dql ^ Rsch), we have also restricted to 



equatorial observers. This considerably simplifies all cal- 
culations without affecting the complete investigation of 
the most significant physical situation, namely the black 
hole in Sgr A*. In fact, it is natural to assume that the 
equatorial plane of this supermassive black hole coincides 
with the galactic plane. In any case, a full investigation 
stepping beyond these restrictions is in progress. 

The first achievement of this paper has been the ana- 
lytical description of the Kerr caustics. At the first order 
in a we find that they are just shifted along the equato- 
rial plane still remaining point-like, while at second order 
they are resolved into typical diamond-shaped figures. 
We are thus able to calculate the position and the exten- 
sion of the caustics for any order of relativistic images (as 
long as we remain in the perturbative regime) . As stated 
in Sect. IV, the strong deflection limit treatment does 
not cover the primary caustic (fc = 1 in our formulae), 
since this caustic is formed in the weak field regime for 
fiS 3> 1 jl^H^. With no regard to the lensing regime, 
the effect of the angular momentum of the deflector is 
similar, with caustics getting a diamond-shaped struc- 
ture and drifting from the optical axis. Whereas magni- 
fication effects due to the primary caustic are very large, 
its size is very small, so that the creation of additional 
pairs of images (which is the most evident manifestation 
of the presence of a non-negligible spin) is very difficult to 
achieve. The significant extension of relativistic caustics 
strongly enhances the cross-section for additional images 
and puts them in a much better position for testing the 
Kerr nature of the black hole. It has to be remarked that 
our perturbative investigation still leaves open the possi- 
bility that metamorphosis may occur at large values of a. 
Though this has been numerically excluded for the pri- 
mary caustic |26j , it is possible that higher order caustics 
develop more complicated structures in a strongly non- 
pertubative regime. 

The second achievement has been the analytical inver- 
sion of the lens mapping near the caustics, which has 
allowed us to draw fascinating pictures of the relativis- 
tic images generated by a hypothetical source close to a 
relativistic caustic. 

However, the most important result has been the pos- 
sibility of doing concrete analytical estimates of the 
size and luminosity of the relativistic images. The 
LMXBs surrounding Sgr A* provide an ideal population 
of sources, which may eventually bump into a relativis- 
tic caustic and generate appreciable relativistic images. 
This is because they are compact sources with very high 
surface brightness in the X-rays, compared to that of Sgr 
A*. This seems not to be the case for stellar sources 
in the infrared K-band, which have a too small surface 
brightness. Using our formulae for the magnification of 
relativistic images, we are entitled to claim that future 
space missions performing X-ray interferometry with res- 
olutions of the order of 1 jias will see these relativistic 
images with high probability. 
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APPENDIX A: RESOLUTION OF RADIAL 
INTEGRALS 

In this appendix we recall the SDL technique used in 
Ref. (171] to solve radial integrals, applying it to the inte- 
grals that appear in the geodesies equations. We rewrite 
them here for an easier reading 



/i = 2 



dx 



R 



(Al) 



(A2) 



First we change the integration variable from x to z by 
the transformation 



1 - z' 



(A3) 



As a consequence, the integration domain [a:o,oo] be- 
comes [0, 1]. 

Then, each of the integrals /i, I2 can be written in the 
form 



/, - y R,{z)f{z)dz 



1 



(A4) 
(A5) 



where the two functions Ri{z) can be easily read by Eqs. 
(|24|l - (|25|l taking into account the Jacobian of the trans- 
formation HA3|I : 



Ri{z) 



2xn 



(l-z)2 



(A6) 



R f^):^ 2a;o xl + {l-z)^{a^-aJ) 

' (1 - z)2 xl - a;o(l -z) + a^l - z)2 > 

Now we consider the expansion of R{z) in a neighbor- 
hood of z = 0. Since z = means x = xq and xo is a 
root of R{x), we deduce that R{z = 0) = 0. Then the 
expansion of R{z) reads 

R{z)=az + /3z^+o{z^), (A8) 

where the coefficients of the expansion are 

Xo [(a - J)2 + g + 2{a^ - - Q)xo 



+4x1] 



(A9) 



13= Xo [(a -Jf + Q + 3(a2 - - Q)xo 



+lQxl\ 

We use this expansion to define 

1 



^Jaz + j3z 



(AlO) 



(All) 



The radial integrals can be split in two pieces 

li — Ii,D + liM. 



h.D = / R^{0)h{z)dz 



(A12) 
(A13) 



hR = J [Rr{z)f{z) - R^iO)foiz)] dz. (A14) 





The first integral gives the result 



2i?, , V^-t-y^T? 

D = — 7^ log j= . 



(A15) 



a, (3 and i?i(0) are known functions of xq, J, Q and a. 
Now we can use the SDL parameterizations ()20|l. H21|l . 
()22|l for all these quantities, so that they become func- 
tions of ^, S (or equivalently e) and a. Then, in the spirit 
of SDL approximation, we keep the leading order in 5, 
which goes as log 6, and the next-to- leading order which 
is constant in 6. Finally, we expand the obtained expres- 
sion to second order in a. 

As regards the integrals h^R, the integrand function is 
regular in the whole integration domain. Sending S to 
zero, the integrand does not diverge. So, this integral 
contributes to the SDL expansion with another constant 
in S plus higher order terms that we can neglect. It is con- 
venient to make the second order expansion in a before 
the integration, in order to have a sum of easily integrable 
functions. We can then add the result of the integral li^j^ 
to the integral h^o, to reconstruct the full SDL formulae 
for radial integrals: 



Ii = -ai log 5 + bi 
h = -02 log 5 + b2 



(A16) 
(A17) 

The coefficients expanded to second order in a read 
4 16 



^1 a£ H 1= 

3^3 27 ^ 27%/3 



(1 + 3r)a (A18) 



0-2 

62 



hi = ai log 6 {V?> — 

-^«'(5-2V3)(l-a (A19) 

^ + ^«e+^(l + 7eV (A20) 
[a2 + 4(1 + 2a2)] log [6(V3- 1 



-2(1 -I- 2a2) log 48- 



3V3 



a^(3\/3- 2) 



27 



a2[19\/3- 12-^^2(^4 _25V3)]. (A21) 



The separation of Ii into ^ and U^r is necessary 
to isolate the term generating the log (5 into an easier 
integral. 
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APPENDIX B: ANGULAR INTEGRALS 

This appendix is devoted to the resolution of the an- 
gular integrals 



Ji ± / -^d^ 

J Ve 

J2 = ± 



(Bl) 
(B2) 



Introducing the variable fi = cos-d, the two integrals be- 
come 



Ji =± 



Jo =± 




(B3) 
(B4) 



where 



bjQ = 



JQ 



a' - Ji 



2a2 



(B5) 

(B6) 
(B7) 



and we have already replaced J and Q with Jm and Qm, 
coherently with the fact that we only retain terms which 
are logarithmically diverging or constant in S (or equiv- 
alently e). 

8^ has two zeros in fi = ±/j,+ . Then the photon per- 
forms symmetric oscillations of amplitude w.r.t. the 
equatorial plane. It is useful to write the explicit ex- 
pressions of /i+ and /i_ in terms of a and ^, using Eqs. 
(fT ^ -l(0 )l and expanding to second order in a. We find 



4ae 8a^£,^ 
^3Vl" 27 
a(4-8e2) 



2a ' 3%/3 

4a2^(12 - 17^2) 
27 ■ 



(B8) 



(B9) 



The oscillation amplitude is y'l — £^ plus corrections 
due to the black hole spin. This is coherent with the 
fact that for a photon reaching the observer from the 
equatorial plane — ±1) the amplitude of the oscillation 
goes to zero. On the other hand, a photon moving on 
a polar orbit — 0) performs oscillations of maximal 
amplitude, touching the poles of the black hole. Now, 
to perform the angular integrals, it is wise to expand the 
integrands to second order in a and then integrate. Then, 



the primitive functions read 



2 . fj. 

arcsm 



3V3 
81^3 



arcsm 



(33^2 - 9) arcsin 



+ 



(1 -^2 _ ^2)3/2 



(BIO) 



3%/3e 
8 

4a2 
81\/3 



arctan ■ 



arctan 



12 - 21^2 



4M(l-a(3M^- 3 + 4^2) 

^2(1 _ ^2)3/2 



' arctan - 



arctan 



Vi - - e 



(Bll) 



The integration limits are the values of fi at the ob- 
server and source position. The observer is at /Xq — 0, 
since it lies on the equatorial plane, while the source is 
in Us = COSIES- Notice that the choice of an equato- 
rial observer leads to a considerable simplification, since 
Fji (0) — Fj^ (0) = 0. Moreover, we have to consider that 
during the photon motion, /i may perform several oscil- 
lations between — and depending on how many 
loops the photon makes around the black hole before es- 
caping. So, we have to add an arbitrary integer number 
TO of integrals covering the whole domain [— /i+ , . The 
final results read 



2TO7r 



3V3 
2a2(12 



4ae ^ a^(llg=^-3) 
2TO7r 



1 



3V3C 
21^2 + ^3 



9 

4a 
3V3C 
20^"*) 



27^2 



(B12) 



(B13) 



where the minus signs hold if the photon initially in- 
creases its latitude and the plus signs hold if the latitude 
decreases initially. 
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